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Abstract 



The phase structure associated with the chiral symmetry is thoroughly in- 
vestigated in de Sitter spacetime in the supersymmetric Nambu-Jona-Lasinio 
model with supersymmetry breaking terms. The argument is given in the 
three and four space-time dimensions in the leading order of the expan- 



^ . sion and it is shown that the phase characteristics of the chiral symmetry 

is determined by the curvature of de Sitter spacetime. It is found that the 
symmetry breaking takes place as the first order as well as second order phase 
transition depending on the choice of the coupling constant and the param- 
eter associated with the supersymmetry breaking term. The critical curves 
expressing the phase boundary are obtained. We also discuss the model in 
the context of the chaotic inflation scenario where topological defects (cosmic 
strings) develop during the inflation. 
04.62.-Fv, ll.30.Na, ll.30.Rd, 98.89.Cq 
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I. INTRODUCTION 



In the scenario of the early universe it is understood that the grand unified theory 
phase is broken down to the phase of quantum chromodynamics and electroweak theory 
through the Higgs mechanism. While Higgs fields are normally regarded as elementary 
fields, it is of interest to consider a possibility that the Higgs fields may be composed of some 
fundamental fermions as in the technicolor model and to see the consequence of this idea 
in the scenario of the early universe. On the other hand the supersymmetry is supposed to 
be a vital nature possessed by the fundamental unified theory and hence the incorporation 
of the supersymmetry in composite Higgs models is principally important. Under these 
circumstances it is natural for us to consider a supersymmetric composite Higgs model in 
the early stage of the universe and to see whether any remarkable effects are drawn during 
the infiation era. 

The Nambu-Jona-Lasinio (NJL) model is a useful prototype model to investigate the 
mechanism of the dynamical symmetry breaking In many composite Higgs models the 
NJL-type Lagrangian is employed to realize the dynamical Higgs mechanism. From the 
standpoint of exploring the unified field theory of elementary particles it may be of interest 
to investigate a possibility of a supersymmetric version of the NJL model. Unfortunately, 
however, in the supersymmetric version of the NJL model the chiral symmetry is strongly 
protected to keep the boson-fermion symmetry and hence the dynamical chiral symmetry 
breaking does not take place . If a soft supersymmetry breaking term is added to the su- 
persymmetric NJL Lagrangian, the dynamical breakdown of the chiral symmetry is brought 
about for sufficiently large supersymmetry breaking parameter A 0. The reason for this 
is simple: The large A implies the large effective mass of the scalar components of the 
superfields so that quantum effects due to the scalar components get suppressed compared 
with that of the spinor components. Thus the model becomes closer to the original NJL 
model which allows the dynamical fermion mass generation. Stating the same substance in 
a different way we realize that the boson by acquiring its mass term forces supersymmetry 
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to make balance so that the fermion mass is generated dynamically. 

The supersjTumetric NJL model with soft supersymmetry breaking is useful to study the 
mechanism of the dynamical chiral symmetry breaking within the framework of supersym- 
metry. If we take the model seriously as a prototype of the unified field theory, it is natural 
to extend the argument to take into account circumstances of the finite temperature 
and spacetime-curvature as in the early universe 

At the inflation era the quantum effect of the gravitation is of minor importance while 
the external gravitational field is non- negligible. Hence we are naturally led to the super- 
symmetric NJL model in curved space. Dealing with the composite Higgs fields is essentially 
nonperturbative and does not accept approximate treatments. Accordingly we try to solve 
the problem rigorously working in a specific space-time, the de Sitter space, which possesses 
a maximal symmetry. The de Sitter space is suitable for describing the inflationary universe. 
As a nonperturbative method we rely on the expansion technique. 

The four-fermion interaction model (which is the basis of the NJL model) in de Sitter 
space has been discussed by several authors and is found to reveal the restoration 

of the broken chiral symmetry for increasing curvature as a second order phase transition. 
The supersymmetric version of the NJL model in curved space was considered by I. L. 
Buchbinder, T. Inagaki and S. D. Odintsov in the weak curvature limit. They found that 
the chiral symmetry is broken as the curvature increases. Their result is in contrast with 
the result in the nonsupersymmetric NJL model. On the other hand the supersymmetric 
NJL model in the fiat space-time has been investigated by several authors in the context 
of dynamical chiral symmetry breaking. 

In the present paper we investigate the chiral symmetry breaking phenomena in the 
supersymmetric NJL model in de Sitter spacetime induced by the varying curvature. The 
situation is considered to be suitable to simulate the phase transition during the inflation- 
ary period. In the inflationary period the universe expands rapidly with increasing speed. 
This phenomenon is often called the de Sitter expansion. Many investigations on quan- 
tum phenomena in de Sitter spacetime have been performed motivated by the inflationary 
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paradigm. The present authors have recently investigated the chiral symmetry breaking in 
the super symmetric NJL model in de Sitter spacetime in ref. 0. In ref. |^ we have mainly 
worked in the case of three spacetime dimensions. In the present paper we extend the pre- 
vious investigation to the model of four spacetime dimensions, and discuss the cosmological 
consequence of the symmetry breaking phenomenon. 

This paper is organized as follows: In section 2 we describe our model of the supersym- 
metric NJL model in de Sitter spacetime. Here formulae to obtain effective potentials in 
de Sitter spacetime in the 1/N expansion method is described. In section 3 we consider 
the case of three spacetime dimensions D = 3. The results in section 3 partially overlap 
the ones in ref. 0. We review the phase structure of the chiral symmetry because the case 
D = 3 is very instructive for considering the case D = 4. The case D = 4 is investigated 
in section 4. In section 5 we consider our model in the context of the chaotic inflation and 
discuss the possible formation of cosmic strings during the inflation. Section 6 is devoted to 
the summary and discussions. In this paper we use the units h = ks = c = 1, and adopt 
the convention (— + + +) for metric and curvature tensors 



II. SUPERSYMMETRIC NAMBU-JONA-LASINIO MODEL IN DE SITTER 

SPACE 

A. Lagrangian for supersymmetric NJL 

In this section we first summarize the basic ingredients of the supersymmetric NJL model 
in curved spacetime and then we evaluate the effective potential in de Sitter spacetime. We 
consider the following Lagrangian for the supersymmetric NJL model expressed in terms of 
the component fields of superfields 0, 

-^(;7 + p)^-^p^ (2.1) 



4 



where is the number of components of boson field and fermion field ip, respectively, A is 
the four-fermion coupling constant, and p is the auxiliary field. Note that in Eg. ( |2.1D only 
relevant terms in the leading order of the expansion are exhibited. 

Since the supersymmetry of the model protects the chiral symmetry from breaking , 
we introduce the following supersymmetry breaking terms, 

6C = -A20t0 _ A^^^t^c _ ^^/^0t0 _ ^2i?0"V, (2.2) 

where R is the spacetime curvature, and ^i, ^2, and A are coupling parameters, respec- 
tively. 



B. Calculation of effective potential 

In this subsection we describe the effective potential in curved space. Our strategy is 
based on the expansion and we obtain the effective potential in the leading order of the 
expansion. The partition function for the Lagrangian C + 6C in D spacetime dimension is 
given by 



+ p)^ - —p'} 



(2.3) 



up to the normalization constant, where pi and p2 are defined by 



p] = + + i,R, 



(2.4) 



with i = 1, 2. Integrating over 0, (p'^, ip, and ip, we find 



Z = J Vp[det{WV, - plT^'ldetiV^V, - pi)] 

r f N 
x[det{;^ + p)f exp[t J d''xV^{-—p^) 



2\i-N 



-N In det{V^V^ - _ - ^^R) 
-N In det(V^V^ - - - ^2^?) 
+N\ndet{/7 + p) . 



The effective action for Large N is written as 



nP]=Seff[p]+0{-), 



(2.5) 



(2.6) 



and the effective potential V{p) in the leading order of the expansion can be explicitly 
calculated such that 

Seff[p = const] 



V{p) 



1 

-i In 



+i In 

+i In 
1 



-( P \ 



det{P + p) 

det{p) 
detiV^'V^ - p^ - - iiR) 

rfet(V^V^ - A2 - ^^R) 
rfet(V^V^-p2- A2-^2i?)- 



rfet(V^V^ - A2 - ^^R) 



-( P \ 



-i{Trln(^ + p) - Tr ln(;e)} 
+i{Tr ln(n - A^ -^li?) 

-Tr ln(n - A^ -^i?)} 
+i{Trln(n-p2- A^-^ai?) 

-Tr ln(n - A^ -^2^?)} 



(2.7) 



where f2 is the spacetime volume and □ denotes V^V^. The effective potential ( ^.71 ) is 
normalized so that V^(0) = 0. Eq.( p^.7D may be rewritten in the form, 



\^ ds tr(x|(/7 + sY^\x) 
Jo 

+i dt {x\{n -t- A'^ - ^iR)-^\x) 
Jo 

+i dt (x\(n-t- A"^ -^2Ry^\x). 
Jo 



(2.8) 



It is important to note that the last three terms on the right-hand side of Eq.( p.7|) are 
related to the massive boson and fermion propagators S{x, y; s) and G{x, y; t). If we define 
the functions, S{x,y; s) and G{x,y]t), which satisfy 

(y + s)5(x,y;s) = ^i^iM, (2.9) 
(n-t)Gix,y;t) = ^-^!^lh^, (2.10) 



-9[x) 

respectively, it is easy to find the relationship : 



Six,y;s) = {x\{;^ + s)-'\y), (2.11) 
Gix,y;t) = {x\{a-t)-'\y), (2.12) 



by making use of the following equality 



d^x^-g{x) \x){x\ = i, (2.13) 
(.1.) ^ (2.14) 



Thus the calculation of the effective potential in the leading order of the 1/A^ expansion 
reduces to the evaluation of the propagators in the de Sitter spacetime, 

V{p) = — +i dsTTS{x,x; s) 
IX Jo 

+i / dtG{x, x;t + A'^ + ^iR) 
Jo 

+i dtG{x,x;t + A"^ + ^2R) ■ (2.15) 
Jo 

The next task is to write down the fermion and the boson propagators in de Sitter spacetime. 

C. Boson and fermion propagators in de Sitter spacetime 

In this subsection we briefly review the boson and fermion propagators in de Sitter space- 
time. The de Sitter spacetime is defined as the maximally symmetric curved spacetime. The 
quantum field theory in de Sitter spacetime has been studied extensively and the propagator 



for scalar fields is well-known |TJ-|T8[. The scalar propagator with mass squared t in the D 
dimensional de Sitter spacetime is explicitly written as 

G{x,y;t) = -i ^^^-^^/^ — :^j^:j^2Fi{a,,bs,c,; 1 - z), (2.16) 

where r(a) is the gamma faction, 2-^1(0, b,c;l — z) is the hypergeometric function, and 

(2.17) 







as = 




bs = 










D 


Cs = 


2"' 



(2.18) 
(2.19) 

with 

z=^, a^=^-{x-yf. (2.20) 

Note that a is in proportion to the geodesic distance between the two points x and and 
r is the radius of de Sitter space, which is related to the Hubble parameter by if = 1/r. It 
should be emphasized that the non-minimal coupling term gives effective mass D{D — l)C,/r'^, 
where R = D{D - l)/r'^. 



The fermion propagator is written as |19,20 



Six, y- s) = {A{x, y- s) + B{x, y- s)a,^r)U, (2.21) 

where U is the matrices composed of the Dirac matrices and 

^ . sr'-^ T{af)T{a}) 
^^'^'^^~\47r)^/2 r(c/) 

x^/l- Z2Fi{af,a},Cf;l- z), (2.22) 

where 

D D , , 

af = j + ^sr, cf = - + l. (2.23) 

We do not present the explicit expression of the invariant function B{x,y;s) because it is 
irrelevant to our purpose of calculating the effective potential. In fact we find 
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TtS{x, x; s) = lim A{x, y; s)TtU 



2^0 



Tt[1] lim A{x,y;s), (2.24) 



with the normahzation hm^^oTrf/ = Tr[l]. 

Now we are ready to calculate the effective potential ( |2.15D . It is instructive to investigate 



the case = 3 because the effective potential can be written in terms of elementary 
functions. We first review the phase structure of the case D = 3 in the next section, and 
then the case D = A is studied in section 4. 



III. ANALYSIS IN 3 DIMENSIONS 
A. Gap Equation 

Following the previous argument, the expression for the effective potential is given by 



V{p) = -^ + i[ dsTTS{x,x; s) 
IX Jo 



+2i / dtG(x,x;t) , (3.1) 

where we set = ^2 = ^ for simplicity. 

In order to write down the effective potential we consider the coincidence limit of the 
propagator, limy^xG{x,y;t). In the coincidence limit, the propagator diverges in general. 
Then we adopt the point splitting method and consider G{x,y;t) around z = 0. For that 
purpose we use the mathematical formula (e.g., pT|]): 

r(a)r(6) ^ , , 

r(c) 

r(a + 6 - c)z''"''~^2Fi{c -a,c-b,c-a-b + l;z) 

+ V! ^2i^i b,a + b-c+l;z) . 3.2 

1 (c — a)L [c — b) 

With the use of this formula the boson propagator is written as 

G{x, y; t) = --^ + -^i/cothTTi/ + 0{z^l'') , (3.3) 



where we defined v = y/tr"^ — 1. Analytic continuation is needed for tr"^ < 1. On tlie otlier 
liand for 2; Eq. (|2.22| ) reduces to 

is is /I r, ^\ tanhvrsr 

y; s) = 5 — - ^ 7 + 

Svrrz^/^ Airr \4 ) sr 

+0{z^'^) . (3.4) 
Inserting ( |3.3| ) and ( |3.4| ) into then the effective potential turns out to be, 

y(p) = -^H / dss{- + s^r^] 

dtucothnu, (3-5) 



27rr Ja^+^r 

where we have set Tr[l] = 4.Q Note that the divergent terms cancel each other, which would 

originate from the supersymmetry. 

The gap equation, dV/di^p^) = 0, reads 

nr /I 9 r,\tanh7rpr 

— + - + ^ 

A \A J pr 

-^p2r^ + 6C - lcoth7ry^p2r2 + 6C - 1 = , (3.6) 

where we defined 

C = e + ^AV. (3.7) 
o 

Note that the solution of this gap equation is completely specified by two parameters ( and 
vrr/A. This means that the phase structure derived from the effective potential is completely 
specified in the plane of these two parameters. 

B. Phase Structure 

We present the typical shape of our effective potential in the case D = 3 in Figure 1 (see 
also [0). In the region labeled by 'S', the effective potential behaves so that the symmetry 



For D = 3 we adopt the reducible representation of the Clifford algebra of Dirac matrices in 
order to guarantee the existence of 75 . 
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is unbroken. In the region labeled by 'B', the symmetry is broken through the second 
order phase transtion. Finally in the region labeled by 'F', the symmetry is broken through 
the first order phase transtion. The boundaries in Fig. 1, which separate the phases, are 
obtained by direct observation of numetical analysis of the effective potential. Analytically 
the condition is described as follows. The solid line in Fig. 1 is found by solving 

dV 



which is explicitely written as 



= 0, (3.8) 

p=0 



y + ^ - v^ec-icothTrv^ec-i = o. (3.9) 

On the other hand, the condition for the dashed line is rather complicated. The condition 
is V{p^) = 0, where p^: = max{pi,p2} when the equation dV/d{p^) = has two different 
solutions pi and p2. 

The branching point in Fig. 1 is of special interest. It is a critical point which divides the 
broken phase into type 'F' and type 'B'. At the branching point C the following conditions 
are found to be satisfied simultaneously: 



dV 



9(p2 



'''' 



p=0 



9(p2)2 



= 0. (3.10) 

p=0 



The conditions are explicitly given respectively by Eq. ( |3.9|) and 

TT^ 1 coth7r-v/6C — 1 

12 2 sinh^ tt V6C - 1 27rV6C - 1 ' 

From Eq. ( |3.11|) we find that ( = 0.290138(=: C*) at the branching point. By substituting 
this value for ( in Eq. (|3.9| ) we obtain nr/X = 0.08306 at the branching point. 

Since we have seen the phase structure of our effective potential, we now discuss the time- 
evolution of the chiral symmetry assuming an inflation background. This consideration will 
make sense if the Hubble parameter (radius of de Sitter spacetime) changes slowly and the 
effective potential in the inflation phase is well described by that in de Sitter spacetime found 
above. As we will see in section 5, the Hubble parameter slowly decreases as the universe 
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evolves during the chaotic inflation in actual. And the curvature of de Sitter spacetime 
slowly decreases (radius r increases). 

First let us consider the case A = 0. In this case ( is equal to ^. If the coupling constant 
A and ^ are fixed, we move from left to right by increasing r in Fig. 1. This indicates that 
the chiral symmetry will be ultimately restored for sufficiently large r. The initial situation 
depends on the value of ^. By direct observation of the gap equation one can easily show 
that if the parameter ^ is kept below 1/4 the effective potential stays in the region 'S' in 
the limit of 7rr/A = 0. Thus the phase transition does not occur for ^ < 1/4. On the other 
hand the effective potential is in a broken phase of type 'F' or 'B' if we take the value of ^ 
above 1/4. Then the phase restoration occurs as the curvature decreases. 

Second we consider the case ^ = and A 7^ 0. In this case ( — A^r^/6. If the couphng 
constant A and ^ are fixed, then the trajectory for varying r on the (C,7rr/A) plane will be 
the one as shown in Fig. 2. By increasing the radius r, we move from the bottom to the 
top along the straight line. For each hne, the value of AA/tt is fixed (a) AA/tt = 100, (b) 
AA/tt = 15.9, and (c) AA/tt = 2. Along the curve for AA/tt = 100 a first order phase 
transition occurs, however, a second order transition occurs for AA/tt = 2. The critical case 
is AA/tt = 15.9. For AA/tt < 1 the phase transition does not occur. 

Thus in the model A 7^ and ^ = 0, the phase of the symmetry is in the symmetric 
phase first, and it is broken as the curvature of the universe decreases as long as AA/tt > 1. 
On the other hand the symmetry is ultimately unbroken in the case A = 0. This case is not 
relevant to our conventional idea of symmetry breaking. Thus we focus our attention to the 
case A ^ and ^ = in the following arguments. 

IV. ANALYSIS IN 4 DIMENSIONS 
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A. Gap Equation 



In this section we consider the case D = 4. From Eq. (|2.15| ) or the gap equation, 
dV/d{p'^) = 0, reads, 

^ + ^TtS{x, X- p) + 2iG{x, X- p2 + A2) = 0. (4.1) 

Here we consider the case A 7^ and C,i = C,2 = 0, because this choice of parameter is 
relevant to the conventional idea of symmetry breaking phenomena as described in section 
3. 

The boson propagators (|2.16|) in the case = 4 reduces to |: 

-iH^ T {3 /2 + iy)T {3/2 - u) 



G{x,y;t) 



(47r)2 r(2) 
X2Fi(3/2 + z/,3/2-z/,2;l-z), (4.2) 



where we defined 



19 t A2 



(4.3) 



While the function A of Eq.(g]22D is 

isH^ r(2 + is/H)T{2 - is/H) 



A{x,y; s) 



(4vr)2 r(3) 



X VT^2Fi (2 + is/H, 2 - is/H, 3;l-z). (4.4) 

Since the propagators diverges in the coincidence limit, 2; 0, we adopt the point 
splitting regularization by keeping z small but finite. In this case it is useful to use the 
mathematical formula (e.g., pl| ) 

JMFi^F,ia,b,a + b-l;l-z) = 
T{a + 0—1) 

l + ia-m-1) t^^{Hz) + ^{a + n) 

+V'(6 + Ti) - V'(^ + 1) -^("' + 2)}^'", (4.5) 



^In this section we use the Hubble parameter H{= 1/r) instead of the radius r to describe the 
curvature of de Sitter spacetime. 
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where (a)„ = T{a + n)/T{a) and ipi^) is the polygamma function. With the use of this 
formula the gap equation reduces to 

{2tx/HY 1 - y/T^ 

A ^ i 

{2 + ip/H)n{2-ip/H)^ 



X {ln(z) + V(2 + ip/i/ + n) + i){2 - ip/H + n) 

-?/'(n + 1) - V(n + 2)}z" 

(3/2 + z/)„(3/2-z/) 



x{ln(^) + V(3/2 + z/ + n) + ^(3/2 - + 

-^(n + 1) + 2) }z" = 0. (4.6) 

The divergent terms of order 0{z~^), which appears in the boson and fermion propagators, 
cancel each other in the effective potential. This would be traced back to the supersymmetry 
of the theory. However the logalithmic divergent terms 0{\n.z) arise in the case D = A. 
Thus our theory depends on the cutoff parameter. Instead of z we introduce the momentum 
cutoff parameter A. We can show that the geodesic distance a = {x — yY /2 is related to 
the momentum cutoff parameter A such that a = A~^. Then we have the relation 

^ = ^- <^-^' 

where we used Eq.( p.20| ). 



B. Phase structure 

Next we consider the phase structure of the model by solving the gap equation ( [4.61) . 
Because of the complexity of equation, an analytic approach would not be useful in the 
case D = 4. Then we adopt numerical method by calculating equation ( [4.61) . As the gap 
equation depends on the cutoff parameter, we do not find a simple scaling relation in gap 
equation ( [4.6[ ), on the contrast to the case D = 3. Therefore we fix the cutoff parameter 
A = lO^^GeV in the following arguments, for simplicity. 
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Fig. 3 shows the phase structure of the effective potential on the (if/A, l/AA^)-plane. 
The curves in the figure show the phase boundary between a symmetric phase and a broken 
phase for A = lO^^GeV, 3 x lO^^GeV, lO^^GeV, respectively. In the right and upper region 
of each curve in Fig. 3, the effective potential behaves so that the symmetry is unbroken. 
On the other hand in the left and lower region of each curve the symmetry is broken. We 
discuss the type of the broken phase, i.e., first order or second order, in the later. 

By decreasing the curvature (decreasing H) with the coupling constants. A, A, and A 
fixed, we move from right to left in Fig. 3. This indicates that the effective potential is in the 
symmetric phase initially, however, the symmetry is broken as H decreases if the coupling 
constant A is larger then a critical value. Let us find the critical value of the coupling 
constant, which depends on the supersymmetry breaking mass A and the cutoff parameter 



A. Taking the limit of Eg. ( [4.6|) , we find the following critical coupling constant Acr, 



1 _ _ A2 ^ 1 
K.~~'{2nYt,n\{n + l)\ 

A2 .J/A2\- 



X 



In — - -^(n + 1) - V^(n + 2) 



(4.8) 



V2A2, 

We show the critical value A^ as a function of A/A in Fig. 4. We see the expected result in the 
figure that the critical coupling constant Acr must becomes larger as the the supersymmetry 
breaking mass A becomes smaller. 

Finally in this section we discuss the effective potential of type of the first order phase 
transition. Similar to the case D = 3 we find that the effective potential behaves so that 
the phase transtion occurs through first order transition for a specific range of coupling 
constants. The region is small and we can not recognize it in Fig. 3. Fig. 5 shows the region 
in Fig.3 under magnification for demonstraition, where we adopted the case A = lO^^GeV. 
In the narrow region labeled by 'F' the effective potential has the shape of the first order 
phase transition. For other cases, A = 3 x lO^^GeV and A = lO^^GeV in Fig. 3, we find 
narrow regions of 'F' in a similar way. Thus the model for D = 4 has the similar phase 
structure as the case D = 3. 

Because there appears many parameters in the case D = A and we can not find a good 
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scaling relation on the contrast to the case D = 3. Due to this fact it is difficult to show 
the phase structure in case D = 4 clearly. As will be discussed in the next section the 
change of the Hubble parameter is rather fast in the chaotic inflation model. Then the 
period that the the effective potential has the shape of the ffist order phase transition is 
very short compared with the Hubble expansion time. Then the cosmological importance of 
the ffist order transition is unclear at present. Therefore we only focus on the epoch of the 
(second or ffist) phase transtion in the inflationary universe and discuss other cosmological 
consequences because phase transitions in the universe predicts topological defects in general. 



V. EPOCH OF PHASE TRANSTION IN THE CHAOTIC INFLATIONARY 

MODEL 

In this section we consider a cosmological application of our model in the chaotic inflation 
universe 



. Cosmological phase transition predicts the formation of topological defects in 
general. In our model U{1) chiral symmetry is broken down and the formation of cosmic 
strings would be predicted by the phase transtion. In general topological defects are harmful, 
as examplified by the monopole problem and domain wall problem. However, cosmic strings 
have been studied with the motivation of the cosmic structure formation. The observational 



confrontation of the cosmic string scenario is still at issue p3 |. 

The investigations in the present paper show that the symmetry breaking may occur 
during inflation era, which is triggered by the change of the curvature of spacetime. As an 
application we consider the phase transition in the context of the chaotic inflation model 
and the formation of the cosmic strings during the inflation. If the inflation lasts sufficiently 
long after the phase transition, the defects will be diluted. However, if the phase transition 
occurs at a suitable epoch of inflation, the dilution of the strings is not completed. Similar 
investigations were given on the basis of a simple model of curvature induced phase transition 



23,|25|. Following their consideration we here discuss the epoch of the formation of the 



cosmic strings in our model during the chaotic inflation era. 
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We assume the chaotic inflation background. To be speciflc we assume that the inflation 
is derived by the field (p with the Lagrangian 

= -V^(j)V^(l) - ]^m'<p\ (5.1) 

Here we set m = lO^^GeV [531, which is the constraint from the amphtude of density 



perturbations for a successful infiation scenario. With the background metric, 

ds"^ = -dt^ + a{tfd^\ (5.2) 
the equations of motion under the slow roll approximation take the form, 

3if0 + mV = O, (5.3) 



An 



where H = d/a is the bubble parameter, and the dot denotes (cosmic time) t differentiation. 
The solution of inflation is well known 

m=<l>i-ir^it-U), (5.5) 



a{t) = tti exp 



(5.6) 



16n 



m = \i^m (5.7) 



3MpL 



R = 12H'' + 6H = ii^rn^(l)^{t) - 2m 
IGtt 



m'<p\t), (5.8) 



where MpL is the Plank mass, and the subscript 'i' denotes the value at the initial time 
when the inflation started. 

In this chaotic inflation model the horizon crossing (during the inflation) of the fluctu- 
ation of the wave length corresponding to the present horizon size occurs at 0(t) ~ 3MpL. 
The epoch of forming cosmic string are estimated as follows. The phase of p does not flx 
soon after the phase transition because of the quantum fluctuations of the p fleld. Following 
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the previous investigations on curvature induced phase transition p4| , p5| , p2| , the phase of p 
fixes at the time, when the condition is satisfied 

H^ = -\imV{p)\ (5.9) 

which we may regard as the time of the string formation during the infiation. We denote 
the value of (pit) at this time as 02- With the use of equation (|5^ ), 02 is written as 



— ^hm\/(p)". (5.10) 



In the previous work p4| , |25| the formation of the cosmic strings during infiation was 
discussed to avoid observational difficulties of the cosmic string scenario. They argued the 
range of the model parameter (coupling constant) to realize the suitable formation of the 
cosmic strings. In this paper we consider the observational possibility and find the condition 
that the cosmic strings are formed during the infiation without complete dilution of the 
cosmic strings so that the cosmic strings may exist in the universe at present time. The 
condition is simply given by 

^<02^3MpL, (5.11) 
v47r 

where we have assumed that the strings are formed at the rate of one string per horizon 
volume at the formation epoch. 

This condition constrains the coupling constant A in our model. From Eq.( ^.9|) , we have 



T + - 7^v/l^2i"i(2, 2, 3; 1 



7r(l/4-/i2) , , , 

+7^^^ ^2^1 3/2 + p, 3/2 - /i, 2; 1 - z = 0, 

(27r)^ cosnp 



(5.12) 



where p = a/9/4 — A'^/H'^, z = iJ^/2A^, and H is specified as 



*=\/3Mr^- (5.13) 
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Fig. 6 shows the formation time of the cosmic strings on the parameter space, A and A. 
Here we set A = lO^^GeV. The result depends on the choice of cutoff parameter A. However 
the quahtative feature does not depend on the choice of A. The region between the two sohd 
curves in Fig.6 satisfies the condition ( p.lll ). The lower solid curve is from the condition 
that the formation of the strings occurs at 02 = 3MpL. While the upper solid curve comes 
from 02 = MpL/-\/47r. The larger coupling constant A leads to the earlier formation of the 
strings during the inflation. And the larger breaking parameter A also leads to the earlier 
formation of the strings. Fig.6 shows this expected result. Thus if we make a choice of 
parameters in the lower region than the lower solid curve, the strings are diluted out by the 
inflation and they would not be observable. | The dashed line, which is almost overlapped 
with the upper solid line, shows the critical coupling constant for the phase transitions (see 
also Fig. 4). Therefore the strings do not form when we make a choice of parameters in the 
region above the dashed line. 



VI. CONCLUSIONS 

In summary we have shown that the curvature of spacetime might be a trigger of phase 
transition during the inflation in a class of model of the dynamical symmetry breaking. 
To be specific we have investigated the phase structure of the super symmetric NJL-model 
in de Sitter spacetime. By evaluating the effective potential in the leading order of the 
expansion, we have examined the phase structure of the chiral symmetry in three- and 
four-spacetime dimensions. We have investigated how the curvature of de Sitter spacetime 
changes the phase structure, and we have found that the symmetry breaking takes places as 
the first order as well as second order phase transition depending on the coupling constant 
and the parameter of the supersymmetry breaking in both cases of three- and four-spacetime 



^If the reheating temperature was so high that the symmetry restored by the finite temperature 
effect, the cosmic string would be produced again. 
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dimensions. Q 

In the case of three-spacetime dimensions, the gap equation reduces to a simple equation 
written by elementary functions. Divergent terms cancel each other and we do not need 
introduce a cutoff parameter. As shown in section 3, the shape of the effective potential 
can be specified by the two parameters, r/A and (, which characterizes the phase structure. 
This model might not have physical meaning, however, it is very instructive to consider the 
case of four-spacetime dimensions. 

In the case of four-spacetime dimensions, the gap equation does not reduce to a simple 
form, which is written by the hypergeometric functions. We have to introduce a (momentum) 
cutoff parameter to regularize the theory. The analysis is tedious, however, we have found 
that the phase structure is very similar to the case of three-spacetime dimensions. Similar 
to the case of three-spacetime dimensions we find a narrow region in which the effective 
potential has a shape of first order phase transitions. In the open inflation scenario a bubble 



nucleation (first order transition) occurs during the inflation [^]. It may be of interest to 
consider the possibility whether our model works as a successful model for the open universe 
or not. However this first order phase transition would not work successfully, because the 
curvature (Hubble parameter) of de Sitter spacetime changes rather fast during the inflation 
as we have shown in section 5. 

We have briefly discussed a cosmological application in section 5. We consider the phase 
transition in the chaotic inflation background, neglecting the back-reaction from the fields 
which derives the phase transition. The strong coupling constant A leads to the early phase 
transition and the subsequent inflation dilute the strings. The weak coupling constant does 
not lead to the phase transition during the inflation. Thus this model may work as a model 
of forming the cosmic strings during the inflation if we choose the parameters in a suitable 



^In the case of the potential barrier is small enough, the conventional picture of first order phase 
transition would be broken down because of the quantum fluctuation of the field. 
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range p^,|25[ . 
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FIGURES 
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FIG. 1. Phase structure of the model for £) = 3 on the (vrr/A, C) plane. In the regions labeled 
by 'S', 'B', and 'F', the shape of the effective potential is of type of unbroken symmetry, of second 
order phase transtion, and of first order phase transition, respectively. The square denotes the 
branching point. 
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FIG. 3. Phase boundary in the case D = 4 on the {H/K,l/Xh?) plane. Here A = lO^^GeV 
is used. We adopted the supcrsymmetry breaking parameter A = 10^^, 3 x 10^^, lO^^(GeV), 
respectively, for each curve in the figure. The right and upper region labeled by 'S' is in a symmetric 
phase for each curve. While the left and lower region labeled by 'B' is in a broken symmetry. 
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FIG. 5. Phase boundary for the model for D = 4 on the (i?/A, 1/AA^) plane. Here we set 
A = lO^^GeV and A = lO-'^^GeV. In the narrow region between two curves labeled by F, the 
effective potential has a shape of type of first order transition. The square is the branching point 
obtained from numerical calculation. 
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FIG. 6. Constraint on formation epoch of cosmic strings on the (1/A, A) plane. Here 
we adopted the momentum cutoff A = lO^^GeV. The upper solid curve reflects the condition 
02 = (47r)~^/^MpL and the lower solid curve corresponds to 02 = 3MpL. The dashed line shows 
the critical coupling constant Acr- 
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